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Abstract 


Operating conditions which yielded the peak 
thermodynamic efficiency (41 percent) for an ETF- 
size MHD/steam plant were previously identified by 
considering only the active region (the primary 
portion for power production) of an MHD channel. 

In this paper, the efforts of the previous study 
are extenoed to include the investigation of the 
effects of the channel end regions on the overall 
power generation. Considering these effects the 
peak plant thermodynamic efficiency is found to be 
slightly lowered (40.7 percent) and the channel 
operating point for the peak efficiency is shifted 
to the supersonic mode (Mach No., Mq 'v 1.1) 
rather than the previous subsonic operation 
(Me 0.9). The sensitivity of the channel 
performance to the B-field, diffuser recovery co- 
efficient, channel load parameter, Mach number, 
and combustor pressure is also discussed. In 
addition, methods for operating the channel in a 
constant-current mode are investigated. This mode 
of channel operation is highly desirable from the 
standpoint of simplifying the current and voltage 
consolidation for the inverter system. This sim- 
plification could result in significant savings in 
the cost of the equipment. The initial results of 
this study indicate that this simplification is 
possible, even under a strict Hall field con- 
straint, with reasonable plant thermodynamic 
efficiency (40.5 percent). 

I. Introduction 


Previous studies^*^ have identified MHD 
channel configurations and operating conditions 
which result in peak thermodynamic efficiency of 
the combined MHD/steam plant. A critical factor 
in these studies was that the MHD channel was 
cooled by low pressure, low temperature boiler 
feedwater. This results in a lowering of the 
steam plant efficiency with increasing MHD channel 
heat loss since it reduces the amount of_regenera- 
tive boiler feedwater heating that can be done 
within the steam plant itself. In these initial 
studies, only the active (significantly power pro- 
ducing) region of the channel was considered. 

This region was taken to extend from a B-field low 
of 4 tesla at the front of the channel to a 
B-field low of approximately 3.5 tesla at the rear 
of the channel (see Fig. 1). However, even though 
the regions upstream and downstream of the active 
region do not produce significant amounts of 
power, they still produce significant voltage 
which requires voltage isolation. Therefore, 
these regions require MHD-channel-type construc- 
tion and hence cooling with low-pressure, low- 
temperature boiler feedwater with the attendant 
possibility of a further reduction in overall 
plant efficiency. 

The results presented in this paper apply only 
to plants of the size of the ETF conceptual design 
(200 Ml^) and to MHD channels whose design re- 
quirements specify cooling with low-pressure, low- 
temperature boiler feedwater. The sensitivity of 
various channel parameters (maximum B-field, dif- 


fuser recovery coefficient, generator load parame- 
ter, Mach number, and combustor pressure) are 
examined under the constraints of a maximum axial 
electric field (Ex) of 2.5 kV/m, a maximum 
transverse electric field (Ey) of 4 kV/m, and a 
maximum transverse current (jy) of 10 kA/m^. 

In addition, the channel performance is investi- 
gated using constant current output in various 
regions of the channel in order to take advantage 
of possible cost savings for the inverter system 
under these conditions. 


II. Background 

The channel analysis is based on a quasi-one- 
dimensional code.^ However, in this study the 
code has been modified to allow the calculation of' 
the voltage drops through the integral boundary 
layer. This is necessary since the original model 
IS inadequate for the channel end regions. In 
addition to the voltage-drop calculations, the 
channel code utilizes tabulated chemical equilib- 
rium properties which are computed separately.^ 
This allows a considerable reduction of the re- 
quired computing time, especially in the region of 
the boundary layers. The channel code is then 
combined with an overall system calculation^ to 
determine the performance of the MHD/steam plant. 

III. Channel Characteristics 


The operating conditions and electrical stress 
constraints used in this study are given in 
Table 1. In Ref. 1, the optimal channel had an 
active length of 12 m. In this study, the end 
regions of the channel are included by starting 
the calculations at the 0.5 tesla point and con- 
sidering the entire length of the channel to be 
16 m which results in a B-field of 0.5 tesla at 
the exit (Fig. 1). At these B-field levels, the 
Hall electric field at the channel entrance and 
exit IS approximately 10 V/m and 40 V/m, respec- 
tively. At these electric field levels, voltage 
isolation IS no longer a serious problem and the 
construction is that of the nozzle or diffuser 
which are cooled by intermediate or high-pressure 
boiler feedwater which does not seriously affect 
the overall system efficiency. 

IV. Results 


Table 2 lists the results for twelve cases all 
optimized for combustion pressure (Pc) and load 
coefficient (K^pi^). These cases are for the ETF 
conceptual design input power level of 540 MW^, 
an air enrichment level of 30 percent O2 by 
volume, and a length of 16 m, as discussed below 

Base Case 


Case 1 IS the result of adding the end regions 
to the initial ETF channel configuration in which 
only the active region (central 12 m) was con- 
sidered.^ In case 1, Pc and have 

been readjusted to obtain the peak overall system 
efficiency for the total length of 16 m. An over- 



all efficiency of 40.3 percent was obtained which 
compares to the 41 percent obtained previously. 
Axial profiles of Ex, Ey, Jy, K, MHD power 
(Pmhd)* 3nd Hall parameter (e) for this case are 
plotted in Fig. 2. 

Variation of Mach Number and Maximum B-Field 

The effect of varying the inlet Mach number on 
the plant overall thermodynamic efficiency is 
shown in Fig. 3 for peak B-fields of 6 tesla 
(cases 1-4) and 5 tesla (cases 5-8). It is seen 
that for the Bjnax = ^ telsa case a maximum 
efficiency of 40./ percent is obtained at a super- 
sonic Mach number of 1.1. This compares to a 
value of 41 percent at a Mach number of 0.9 ob- 
tained for the original 12-m channel calcula- 
tion. ^ The increase in optimum Mach number for 
the present calculation is primarily the result of 
the fact that a lower average B-field exists in 
this longer channel. As was shown in Ref. 1, the 
lowering of the B-field results in a higher opti- 
mum Mach number. This is again born out as can be 
seen from the results for the B^ax = 5 tesla 
cases 5-8 where the maximum efficiency of approxi- 
mately 40.4 percent occurs at a Mach number of 
1.25. 

In all of the above cases for Bmax = ^ telsa, 
the diffuser pressure recovery coefficient, Cg, 
was held constant at 0.46. A sensitivity study 
was then conducted by considering Cg = 0.4 and 

0. 55. The result of this variation is shown in 
cases 9 and 10 of Table 2. It is seen that the 
sensitivity in the overall efficiency is about 
±0.2 percent over this range of recovery 
coefficients. 

Constant Faraday Current Operation 

ETF conceptual design studies have shown that 
considerable simplification of the consolidation 
and inverter circuitry can be obtained by 
operating the generator at constant Faraday cur- 
rent, ly = Jy Ae, where Ag is the effec- 
tive electrode area. Cases 11 and 12 of Table 2 
show the effect upon performance of attempting to 
operate various regions of the channel in this 
constant current mode. Because of the direct 
relation of Jy to the axial electric field, 

1. e., Jy 'x* oEx/B, the current level is in 

turn limited by the Ex constraint— especially 
at the downstream section of the channel where the 
electrical conductivity (a) drops and the Hall 
parameter (b) is high. The results of case 11 in 
which an attempt was made to operate the entire 
channel at a constant current output of 
Iy,max 2.9 kA/m (where Iy,max ts the 
maximum transverse current/unit length) are shown 
in Fig. 4, The axial variation of ly, Jy, 

E and B are depicted on the figure. Nearly 
one half of the channel can be run at constant 
current but at X 8 m the Hall field limit is 
reached. 

In an attempt to further simplify the inverter 
system, i.e., to eliminate the continuously vary- 
ing current output at the latter half of the chan- 
nel, a two-level case (Iy,max 3.0, 2.5 kA/m; 
case 12) was investigated! The results are shown 
in Fig. 5 and are not significantly different from 
those in Fig. 4. Using two levels of constant 
current does not extend the length of constant 


current appreciably and does not render signifi- 
cant simplification of the current consolidation 
system over case 11. 

Preliminary calculations show that if 
^x,max relaxed to about 3.5 kV/m, almost the 
whole channel can be controlled at one level of 
constant current but the plant thermodynamic effi- 
ciency drops about one point relative to case 11. 
Nevertheless, both cases 11 and 12 show reasonable 
performance for partial constant-current control 
with a plant thermodynamic efficiency (nth) 
about 40,5 percent. 

Concluding Remarks 

This study of the MHD channel for a 
540 MWth, ETF-size plant, 16 m long, including 
the end sections suggests: 

1. The best performance is obtained in the 
supersonic mode, even for B^^^ax = ^ tesla 

(n 40.7 percent. Mg 'v 1.1). 

2. Lowering B^ax approximately 5 tesla 
does not severely lower the performance 

(nth = ^0*4 percent at Mg 'x. 1.25) and could 
result in a reduction of the magnet size of up to 
40 percent. 

3. The overall performance is not too sensi- 
tive to diffuser pressure recovery coefficient. 

For a Cg range of 0.4 - 0.55, the overall 
efficiency changes by ±0.2 percent. 

4. A simplification of the inverter system 
IS possible. Good channel performance 

(nth percent) can be obtained using a 

partial constant-current mode, even under a strict 
electrical constraint on the Hall field 
(^x,max 2.5 kV/m). If this constraint could 
be relaxed to about 3.5 kV/m, the entire channel 
could be controlled by a constant-current opera- 
tion; however, the penalty in plant efficiency 
might not be negligible. 
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TABLE 1. - THE OPERATING CONDITIONS AND ELECTRIC STRESS CONSTRAINTS 


Thermal input, MW 

Coal type 

Moisture content in coal, percent , . . . 
Oxidizer (O 2 enriched air), vol. percent 
Oxidizer preheat temperature, . 

Seed, 1 percent-mass, K 

Total flow rate, kg/sec 

Combustor heat loss, percent 

Combustor slag rejection, percent .... 

Oxygen stoichiometry, percent 

Poly. eff. of cycle compressor, percent . 
Specific power of ASU, kW-hr/eo. ton O 2 • 
Ambient pressure (atm) and temperature (*F 
B-field 

Axial electric field limit, kV/m .... 
Transverse electric field limit, kV/m . . 
Transverse electric current limit, kA/m^ 


540 

. Montana Rosebud 


. . 5 

. 30 
1100 
K 2 CO 3 
. 133 
. . 5 
. 65 
. 90 
89.8 
, 204 


0.89; 42 

National Magnet Laboratory (NML) 
preliminary design (Fig. 1) 


2.5 

4.0 

10 




TABLE 2 - PERFORHANCE OF THE 16-fii CHANNEL FOR ETF WITH END REGIONS 
(Og ^ 30 percent. Bmax S.6 tpsla, Hg ■>. 0 9 1 35, Cp 0 4 - 0 55) 


Case 

no. 

®»ax» 

tesla 

Mach 

no 

Diffuser 

coefficient 

rotnbustor 

pressure, 

atm 

mo 

power, 

HWe 

Channel 
heat loss, 
HWt 

Cycle 

compressor 

power, 

MW 

02 

compressor 

power, 

MW 

Bottom 

cycle 

input, 

MUt 

Steam cycle 
efficiency, 
percent 

Steam cycle 
output, 
MUe 

Cycle coal 
input, 
MWt 

— 

Gross 

ac 

power , 
KUe 

Thermal 

cycle 

efficiency, 

percent 

jHI 

mm 

0.9 

0.46 

4.4 

85 06 

28 48 

21.98 

11 62 

430 47 

38 76 

133 25 

539 81 

217.5 

40.29 



0.95 


4 4 

86 10 

27.79 

21 98 

11 62 

429 13 

38 80 

133 03 

539 21 

218 3 

40 43 



1.07 


4 4 

87 39 

25 42 

21 98 

11 62 

428 14 

38 96 

133.19 

539 81 

219 7 

40 70 


■ 

1.15 


4 3 

86.29 

23 99 

21 60 

11 62 

428 86 

39 05 

134 23 

539 81 

219 7 

40 69 


n 

1 05 

0 55 

4.4 

81.18 

26 18 

21 98 

11 62 

434 35 

38 91 

135 39 

539 81 

215 9 

39.97 



1 15 


4 ? 

8? 95 

24 10 

21 31 

11 62 

431 91 

39 05 

135 79 

539 81 

217 9 

40 37 



I.?5 


4 0 

81 39 

21.67 

20 42 

11,62 

432 57 

39 20 

137 20 

539 81 

216.1 

40 40 



1.35 


3 7 

77 52 

19 56 

19 17 

11 62 

435 20 

39 51 

140 27 

539 81 

217.0 

40 20 


6 


0.40 

4.5 

83.90 

28 96 

22 36 

11 62 

432 10 

38 71 

133 40 

539 81 

216.4 

40 08 

10 



0.55 

4 5 

87.08 

27 90 

22 36 

11 62 

428 82 

39.79 

132 39 

539.81 

218 6 

40 49 

11 




4 2 

86.26 

28 22 

21 21 

11 62 

428 50 

38 77 

133 30 

539 81 

218.7 

40 51 

1? 


■ 


4 2 

86 51 

27 95 

21 21 

11 62 

428 25 

38 79 

133 27 

539 81 

218 9 

40.56 
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AXIAL DISTANCE - X, m 


Figures. - Axial profiles of E^, Iy,Jy, and p for two- 
level constant-current control/case 12. (ly - 
3. 0/2. 5 kA/m; = 2. 5 kV/m. ) 
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